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ABSTRACT 



This thesis presents a survey of techniques which can be 
used to suppress sidelobes in antenna array systems. Methods 
applicable to fixed arrays include direct null positioning, 
element spacing adjustment, perturbation of complex pattern 
function zeroes, and pattern synthesis. Methods applicable 
to signal processing arrays include performance optimization, 
which can also be extended to permit system self-adaptation 
to changing noise conditions. 
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I. 



INTRODUCTION 



This thesis is intended to present a survey of techniques 
for the suppression of sidelobes of antenna array patterns. 

The use of multi-element sensor structures in seismographic , 
acoustic, and electromagnetic systems greatly extends the 
range of control over the resulting signal transmission or 
reception performance. A single element restricts process- 
ing capability to a single signal; a system of N elements 
provides the opportunity to independently manipulate N sets 
of operating parameters in such a way as to improve various 
suitable measures of system performance. 

Initially we present the basic characteristics of arrays 
and define two fundamental performance measures; signal-to- 
noise ratio, SNR, and directive gain, G. Methods of sidelobe 
suppression applicable to fixed structures are then discussed , 
followed by methods which optimize performance criteria in 
processing systems. 

It should be noted that while the methods described 
encompass a large portion of current theory they certainly 
do not exhaust the field. The capabilities of large-scale 
computer-controlled signal processors, the multitude of use- 
ful array geometries, and the variety of suitable performance 
criteria insure that this field will remain a productive area 
of theoretical and practical research. 



4 



II. BASIC ARRAY PRINCIPLES 



An "array" may be generally defined as an arbitrary 
collection of sensors arranged in an arbitrary spatial pat- 
tern used to simultaneously detect the local states of an 
impinging field or superposition of fields. Historically, 
the concept of an array has been nowhere as general as this; 
arrays have evolved from the extension to additional numbers 
of elements of fixed geometrical patterns, such as the circu- 
lar arrays of underwater sensors, or the linear and planar 
arrays of radio antennas. Theories of array behavior have 
shown a corresponding evolution from specific, regular 
algebraic forms, such as the complex polynomials of Schel- 
kunoff [13] or the Tchebyscheff polynomials of Dolph [7] , 
to the stochastically controlled self-adaptive processors 
of Mermoz [11] , Widrow [17] , and Griffiths [10] , where the 
very indefiniteness of observed signals provides the means 
of prescribing the appropriate processing technique. 

A. PATTERN FUNCTIONS 

The elementary N-element linear array. Figure 1, has a 
total far-field radiation vector 



where V is the radiation vector of each element, a is the 
o n 

complex excitation of the n-th element, S(0) is the space 
factor of the array, as defined by Ramo and Whinnery [12], 




N-l 
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Figure 1. N-Element Linear Array 
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and k is the wave number. The radiation intensity is 



K = K 1st 2 , 
o 



Thus, in systems where the principle of superposition is 
applicable, the field patterns of an array can be arranged 
as a product of the element pattern times a function of 



If a n is a constant in both amplitude and phase, the 
pattern has a principal maximum in the X-Z plane in the 
direction given by 

nkdsinO = 0 

0 = 0 , 

and nulls in the directions given by 

NkdsinO = 2nrf, n<N 



Having chosen a specific array geometry and excitation, it 
becomes apparent that the response of the array in directions 
other than the central direction is no longer under arbitrary 
control; an interfering signal may be located between the 
nulls resulting in degradation of the performance of the 
system. 

The concept of the pattern function can be extended to 
two or three dimensions, using a general coordinate system 
as shown in Figure 2. In the case of regular planar arrays, 
the space factor is 



2 

angular direction. Isl is called the array factor. 





m n 
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Figure 2. General Array Geometry 
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where (9,0) are the direction coordinates of beam orientation, 

and d and d are the element spacing in the X and Y dir- 
x y 

ections, respectively. The general 3-dimensional array has 
a space factor 

/+,. _ jkr coset 

S(9,(?) ~ 2i a n e n n 
n 

where coso^ = sin9sin9 n cos (0-0 n ) + cos0cos9 n , ^ e n '^n' r n^ are 
the position coordinates of the n-th element and c* n is the 
angle between the "look direction" (9,0) and the element 
position vector. 



B. SIGNAL- TO- NOISE RATIO AND GAIN 

Signal-to-noise ratio, SNR, is a critical parameter in 
many receiving systems, such as echo-ranging systems and 
most amplitude-modulation systems. Similarly, gain, or 
directive gain, G, is a critical signal transmission param- 
eter. In considering problems of antenna sidelobes, it has 
been suggested that off-axis signals may cause unwanted 
interference. Conversely, adjustment of the antenna pattern 
to lower sidelobes may degrade on- axis performance. Techni- 
ques of optimization to be described will maximize SNR or G, 
thus achieving a balanced compromise between sidelobe level 
and main-beam degradation. In general, SNR is defined as 



SNR 



output signal power without noise 
output noise power without signal 



This parameter is meaningless without appropriate definition 
of the processing system. A general system is shown in 
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Figure 3. The output signal is 

S = Ts (t +t ) * h (t) 
o m m 

m 

where the incident signal is assumed to have the same struct- 
ure as a function of time at all sensors, T is the time 

m 

delay of the signal at the m-th element, h m is the impulse 
response of the m-th element, and * represents the convo- 
lution operation. The signal power is then 



m — CP 

where the computation has been expressed in the frequency 
domain with signal power spectral density W (00) , filter 
transfer function H (CD) , and radian frequency variable CO. 
The noise output is 

n = ^>n(t + X ) * h (t) , 
m 

so the resulting noise power is 



P 

o 



m 



?/■» 



(CO) H (0)) H* (oO) e^°\ 



m 



n 



m 



60>/2Jt. 



— CO 

An ultimate objective is to manipulate the gain and phase 
shift of the filters in a frequency-independent manner so 
the filter response can be represented as 

H m (0)) = K^e^mH (CO) 

with H (CO) the frequency-independent factor. Then, 



P s " 22 K m K n e:i " ^ n) / w s (a)) l H(a;) l 

m n -os' 



m 



V 60/2 Tt 



and 
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Figure 3. General N-Channel Linear Processor 
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f 0 = HVn eUrpm ' |liMl 2eja)(T m - Vacj.y2n 

m n -on 

Finally, SNR may be expressed in matrix form as 

- K^GK 



SNR 



K T YK 



where 



K = 



V 5 ^ 0 

i 4> 

K, e J Tl 



K M-1 S 



J ^M-1 



K f “ Adjoint [K] 

G = t<W' %r>‘ /W S (“)|HM| 2 e J ®(t m -\) sul/2n 
Y = ty mn i- y ran = /w o (0» IhmI^X - ' r n ) a ^ /2 n 

— CO 

The resulting formula, a ratio of two quadratic forms, pro- 
vides unique properties for optimization. 

Directive gain, G, as defined by Silver [15] is 



G 



power density in direction (0 ,(J? ) 
average power density 



p < e o' < fo> 



2n 7t 

j P (0 , 4>) sinGdGd <P 



Since, for arrays, 

P(0,4>) = p o (Q, ( P) |S (0 , 4>) | 2 

where P Q (0,'|)) is the element power pattern and S(Q,fi) is the 
space factor of the array, the gain becomes 
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p < e o'V 



, 2JT TC ~ 

hx f f p o (0 '^ l s ( 0 '^r sin e d e d ^ 

•X) J o 



Thg array factor is 
2 _ 

n m 



jk(r cos oi -r cos oc ) 



S(0,4»| “ =25> n a m e ]K{r n C ° Sa n 



m 



m 



where the a^ now represent the weighting factors and linear 
filter of the receiving system. 



a = K e 
n n 



j &n 



H (0)) 



Also, the transmitted power density in direction (9 q ,<P o ) is 

p ^o‘V ‘ P o (e o‘V ls(e o'V 12 

so the gain becomes (Cheng [4] ) 

„ . p o< 9 o^o>|s(e o ,(?> 0 )| 2 

G 27T Jt 

3 ft f J P Q (e,<J>) | S (0 ,0) | 2 sin9d0d4> 



o o 



jk(r cos &C -r coset ) 



P (9 ,0 )]>;>a a e JMr n COSa no 
o n m 



m 



mo 



n m 



2 JT TC 



yf f P Q (8 , P) J> n a m e- * ' < r n eoSlX n- r m cosg m > sinedfld^ 



o o 



n m 



= K^AK 

f 

K BK 



This is, again, a ratio of quadratic forms, with K the vector 
of element weights , and 

2JT Jt 



^ , , _ 1 f f _ jk (r cos ct -r cosCC ) 

B ~ [b ] , b ~ / / P (9,(p)e J n n m m . 

1 mn J mn 431 J J o T 

o o 

sin9d9d<£> 
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A = [a ] = DD 
1 mn J 



D = 



jkr coscC 
e J o oo 



jkr.. cosOf 
i 1 
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ikr,, .coso:,,. 
e J M-l (M-1)0 



C. SIDELOBE SUPPRESSION PRINCIPLES 

The most elementary technique of interference elimination 
is to position a null in the interfering direction. This 
method is restricted to situations in which all other 
noises are insignificant and where appropriate care can be 
taken in system alignment to allow for the sharpness of 
nulls. Positioning a null means adjusting the relative 
element phases so that ‘complete cancelation occurs in the 
null direction. Since maximum signal reception or radiation 
requires a different phasing arrangement, the technique leads 
to sacrifice in performance. The suppression of wide sectors 
of noise generally requires more than null positioning. Here 
it is desirable to specify some acceptable level of gain 
suppression in the sector. This result can be achieved by 
adjustment of both element gain and phase. An artificial 
method of gain adjustment involves modifying the positions 
of the elements in a regular array structure. The complex 
polynomial technique of Schelkunoff permits both direct 
positioning of zeroes corresponding to desired nulls, and 
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clustering of zeroes to reduce sidelobe levels over an 
angular sector. 

Optimization processes permit use of the total array 
environment to adjust the element filters to the optimal 
values. Thus, background noise, sectoral noise, and narrow- 
beam interference may all be considered. In addition, the 
observation of the noise environment at the sensors can 
provide information which can be used in a real-time sense 
to readjust the filters, thus permitting the system to con- 
tinually adapt itself to the interference environment. 
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III. FIXED ARRAY SYSTEMS 



Fixed array systems are those in which the parameters 
of the array, such as element position, gain, and phase, are 
determined and the array is constructed according to the 
resulting specifications. Methods of sidelobe cancelation 
in this situation include direct null positioning, spacing 
adjustment, perturbation of complex pattern function zeroes, 
and pattern synthesis. 

A. DIRECT NULL POSITIONING 

SchelKunoff [13] showed that pattern functions for linear 
arrays can be represented by complex polynomials on the unit 
circle 

f(Z) -Zaz", z =e0Msine - <t» 

m=0 m 

where f is a constant phase shift per element and a^ is the 
complex weighting coefficient for the m-th element. Zeroes 
of f (Z) lie on the unit circle and define the positions of 
nulls of the pattern. An M-element array will therefore 
have M-l roots of f (Z) , determined by the a , and conversely, 
the a m can be selected to independently position M-l nulls. 
The resulting design will require individual element feeds 
capable of both amplitude and phase control. 

The concept of virtual couplets provides a null position- 
ing system described by Davies [6] . An M-element array has 
a pattern which can be interpreted by the pattern 
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multiplication principle as the product of a couplet, or 
two-element pattern, and an M-l element space factor. This 
reasoning can be extended to its ultimate reduction of the 



array to a collection of couplets. The couplets must have 
the same relative phasing at each stage. The resulting 
structure for a four-element array is shown in Figure 4. 

In this case, there are three independently-controlled phase 
shifts, which permit control of three pattern nulls. 

B. SPACING ADJUSTMENT 

This method, due to Strait and Cheng [16] uses adjust- 
ments in the spacing of the linear array elements to produce 
results similar to amplitude weighting of the previous 
section. A particular advantage, however, is that element 
gains can now be fixed at values satisfying other require- 
ments, such as fixed-power matched drivers. For a linear 
array of N elements with N even and cophasal excitation, 
the space factor is proportional to 




For a variation in spacing v, 




n=l 




kdsin© + 



f a 2 i ) kVSin ®] 
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Figure 4. Null-Positioning Circuit 
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where kdsinQ, 




For c^l^ sufficiently small, 

Sv -2x n (= os fr Li )T - c n ^K sin ^( 2a i i )]- 

Since sin^ ^\f/\ 

S v «Zl n - -jf- sinjsin 

«2l n [cos(22nijf + -jf (cosptijf- cospj^f}_ 

It is then apparent that, given a desired space factor as a 
cosine series with arbitrary element weights, the series S 
can be equated term by term and the factors c n computed in 
such a way as to hold I n constant. Such a result would 
permit the design of transmitters for optimum, constant 
operating power characteristics, and would also avoid the 
complications of individual element amplification and phase 
weighting. 

C. PERTURBATION OF ZEROES 

This approach involves making small adjustments in the 
locations of zeroes of the complex array polynomial. A 
similar technique has been used by Baker [1] , to modify 
features of another otherwise suitable pattern. Since the 
space factor is 
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N-l 

s(0) = TT 

n=l 



(Z-Z ) 
' n 



the array factor is 



P(e) 



S(9) 



2 _ 



N-l 

TT 

n=l 



-a 






2e n |cosh a^ 



n 




where Z - e^ , Y' ~ kd (sinQ) - <p, 

and Z n = e a n are the roots of the polynomial. Since 

it is assumed that a is not necessarily zero, the root at 
Z^ is not necessarily a null of the pattern. The basic idea 
is to start with a known, approximately satisfactory pattern, 
P o (9), and to construct a desired modification of the pat- 
tern, as shown in Figure 5. Variations may now be made in 

a and b to minimize the excess amplitude of the sidelobes 
n n 

in the sector of interest, or to minimize the mean-squared- 
difference between the -desired pattern and the actual pat- 
tern. Normally, it would be necessary to define the refer- 
ence pattern over the entire range from 9 = 0 to 9 = JT. 
Otherwise, adjustments resulting from differences in a given 
sector may cause reduction of performance in the main beam. 



D. PATTERN SYNTHESIS 

This procedure is also based on the theory of Schelkunoff. 
A desired space factor may be defined by an arbitrary 
function, F(9), over the range 9 = 0 to 9 = Tt, or its 
equivalent function of f (Y) , where 

kdsin9 - $ 
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P(e) 




SECTOR TO BE MODIFIED 



Figure 5. Pattern Modification 
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If the array has an odd number of elements and the coeffici- 
ents equidistant from the ends of the polynomial are required 
to be conjugates, the space factor becomes 

N 

S(Y) = X (A cos ftf + B sin n^|/) 
n = 0 

where the unknown and B n are related to the polynomial 
coefficients. Then, the desired space factor may be expanded 
as a Fourier series so that 

f (\J/) = X (P cos n Y + sin n ^) 

n=0 n n 

and the coefficients of the two series can be equated. 
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IV. LINEAR PROCESSING ARRAYS 



Two types of signal processing systems utilize individual 
filtering, or weighting, of each array element to optimize a 
performance characteristic of the array. In one method the 
signal to noise ratio, SNR, is maximized. Two general 
approaches to SNR maximization have been considered by 
Mermoz [11] using N-channel matched filters, and by Cheng 
[3,4] using the properties of Hermitian matrices. In the 
second method the directive gain of the array is optimized, 
subject to constraints which may include specification of 
the desired location of pattern nulls. The objective of SNR 
maximization is to obtain a set of filters which result in 
maximizing knowledge of the presence of a signal when the 
signal and noise are additive. The maximization techniques 
require knowledge of the structure of noise in the surrounding 
environment; for the matched filter approach, the noise must 
be observable by measurement at the array elements in order 
to determine its statistical properties. This has the advan- 
tage of permitting adjustment of the filters on the basis of 
the noise measurement if the noise environment changes; 
however, it also requires a capability to process direct 
noise information separately from filtered signal plus noise. 
The matrix approach is equivalent to the matched filter 
approach in receiving systems; it permits constraints on null 
locations to be included in the optimization. 
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A. SNR OPTIMIZATION 

It was shown that the signal-to-noise ratio can be 
represented by 

SNR 

KYK 

where K is the vector of element weights and G and Y are 
square matrices relating the relative element time delays 
and the filter, signal, and noise power spectral densities, 
Cheng [3] has used the resulting function to show that if 
G can be represented by 

G = F o P o • 



where 



F = 
o 



-ikr cos CL 
e J o oo 

e- Jkr l cos “o 1 



-ikr w , cos cC . 
e J M-l (M-l) 



then the maximum SNR is 



SNR “ Y 1 F 
max o o 

and the maximizing filter weights are 

K m = Y _1 F . 

M o 

The vector F consists of the relative phase functions of 
o 

the signal at each element of the array. Normalizations 
have been assumed so that the results are valid only within 
a constant factor. Mermoz [11] has used the correlation 
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properties of signals and noises at the various elements to 
arrive at the result that the optimizing filter, H n (G0), is 
given by 

C (CO) H (03) = S *(03) 
mn n m 




where C mn is the noise cross-power spectral density of the 
m and n-th elements and S m is the conjugate of the signal 
spectrum. It has been assumed that the signal input to the 
filters has been delayed to compensate for the relative delay 
of the signals; this operation was previously performed by 
the vector F . The result, then, for optimal performance is 
that 



a) The signal should be delayed by a combination 
of signal time delay and noise correlation; 

b) The signal should be attenuated in proportion 
to the noise correlation; 

c) The result should be linearly filtered by a 
matched filter, H [CO) = S* (00) . 

Combining these results, the complete filter system is given 
by 



H (CO) = (C 1 F ) S * 
m' o m m 



For the M-element system, F q provides the M synchronizing 
delays, C is the MxM matrix of noise correlations (in the 
frequency domain) , and H is the Mxl vector of channel 



filters. S is lxM. If the noise environment is a single 
coherent signal (jamming) from a certain direction, the 
optimization problem reduces to adding the element outputs 
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in such a way as to maintain synchronization of the desired 
signal, while achieving maximum cancelation of the inter- 
fering signal. A two-element monochromatic situation as 
shown in Figure 6 illustrates the application of these 
filtering concepts. For the phase and time-delay reference 
at element 1, the relative signal and noise delays at element 
2 are 

^s = C (d Sin V 

\ = c (d sin V 



If filters and H ^ are assumed to consist of a signal- 
synchronizing delay followed by an inverse noise covariance 
filter, C ^ , then 







i 


|l l| 


~1 0-1 




1 


1 







L J 


C 11 C 12 










e^s . 




*-l -1 




1 


2 








_ C 12 C 22 _ 







In this formulation, amplitude normalization has been 
assumed. Since the noise is coherent and monochromatic. 



c" 1 = c" 1 = 1 

C 11 c 22 1 / 

-l _ *-i , (T - X ) 

C 12 ’ C 21 + eJ « s 



Accordingly, with some algebra. 



" V 

eJ WT s (l-e^N - V) 
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SIGNAL 




Figure 6. 2-Element Filter System 
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For signal and noise inputs A and N at element 1 and Ae 
and Ne-^n at element 2, the filter outputs are 




V 1 = (A+N) Jl - e Ssr q 



V 



2 



(Ae - -*^^ + Ne -^\l) (e-^^s) |l 





Addition of the filter outputs results in the final output. 



V " V + V 0 “ 2A “ 2A cos CO (T - T ) 
o 1 2 Ns 

= 2A [1 - cos 00 {T - Z )] . 

N'S 



Certain comments are appropriate concerning the procedure 
followed before evaluating the result. Under the assumption 
of coherent monochromatic noise, the covariance matrix, C, 
is not invertible, since 

det C = 0. 

This implies that an infinite SNR is obtainable, which means, 
as stated by Mermoz [11] , that the noise is completely can- 
celable. However, the construction of phase-shifting 
filters in the desired structure can still be carried out. 

The procedure also assumes that relative attenuation 
between the sensors is negligible, which would be approxi- 
mately true in practice. The resulting output signal 
contains no noise component; therefore the process of 
cancelation is valid regardless of the angles involved. 
However, the signal amplitude is reduced by delay of the 
noise relative to the signal. If the noise and signal are 
from the same direction, (.^ = , and the output is zero. 

In fact, due to the broadness of the cosine function near 



28 



i - i ~ 0, resolution of the signal and noise will be 
Ns' 

difficult over a large angular range. The resulting output 
signal, V q , should be passed through a matched filter in 
the general case; for monochromatic signals, this would be 
a high-Q resonant filter. 



B. GAIN OPTIMIZATION 

Maximization of directive gain, G, is accomplished 
using the same characteristic of quadratic forms as used 
for SNR. Thus, if 

G = K f AK , 



then 



and 



K f BK 



f -1 

G “ D B D 
max 



K* = B ' 1d ' 



where K m is the vector of maximizing filter weights, A, B, 
and D are as previously defined. Since the elements of B 
are 



2ft ft 



b = 7 Ftt f f p 

mn 4JI J J o 



jk(r cos -r cos oL ) 



n 



n m 



m' sinGdGd^ 



o o 



it is apparent that maximizing G is equivalent to maximizing 
SNR under the assumptions of monochromatic operation with 
isotropic elements in omnidirectional noise. 

It may be desirable to utilize constraints in the opti- 
mization of SNR or G. In the case of transmitting systems, 
radiation in a certain direction may cause excessive inter- 
ference; for receivers, a fixed noise source may be more 
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easily handled by simply requiring a null in its < 

rather than designing the correlation processes o 

needed. Cheng [3] , and Drane and Mcllvenna [9] h; 

lined procedures to accomplish the constrained ma: 

based on transformation of the A and B matrices to an 

abridged form which includes the effect of the constraining 

relations. 

C. ADAPTIVE SYSTEMS 

The signal processing systems previously described can 
be designed for a static mode of operation in which the 
filter weights are predetermined. The array may accordingly 
be adjusted for a fixed main beam direction and fixed noise 
distribution or for scanning beams. The full capabilities 
of the correlation principles used are achieved when the 
design is extended to time-varying situations in which the 
sidelobe structure depends on the orientation of the main 
beam or the noise environment is subject to random variation. 
In such cases, the correlation measurements can be made on a 
real time basis, thereby keeping the system operation in 
nearly optimal condition. The time varying characteristics 
are usually assumed to vary slowly; otherwise, the statisti- 
cal model of the noise may be over-complicated. The main 
problem in developing a self-adaptive arrangement is con- 
struction of the correlation functions. One general scheme 
for control of the filters is shown in Figure 7. Operation 
of this system, as described by Brennan, et. al. [2], is 
defined by 
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ARRAY 




Figure 7. Adaptive System for SNR Maximization 



W = G(F q -Z) 

Z = Q (V*W T V) 



where F q is the steering vector and Q is the Laplace trans- 
form of a linear filter, and W is the vector of weights. 

The desired covariance matrix is 



m ij 



T _ N 

But W V " V " >_ W.V. 

o i i 

i~l 



= E 




so 



Also , 




E 





= E 



MW = 





\i 














N 

z 

i=l 


* 

V, W.V. 
1 11 




< 
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• 

• 
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N 

z 

i=l 


★ 

V x • W.V. 

N 11 






— 




a 



and since 




it is clear that 
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E 




= MW. 



The weights are thus seen to be controlled by the difference 
between the steering vector and a filtered function of the 
product of covariance and weights. In the first-order 
case , 



Q = 



1 

VB+ 1 



and in the second order case 



\\s 2 + + \) s + 1 

Here T represents the eigenvalues of the smoothing process, 
Q, not to be confused with relative signal and noise delays. 
This form of filter control leads to the optimal weights 
which maximize SNR. ' Widrow [17] has considered another 
procedure shown in Figure 8 in which the system adapts to 
match a desired signal by minimizing the mean squared error. 
The pilot signal can be injected internally, or added to the 
external signal and noise field. Doyle's [8] comparison of 
the SNR-maximization and signal-matching systems indicates 
that the pilot-signal technique requires high power to 
achieve adequate "lock-on". 



D. SPATIAL SIGNAL DISTRIBUTIONS 

The systems previously described have modeled signals 
and noise as observable functions of time for which time- 
correlations can be constructed; the external, spatial 
properties of the signals have been assumed to be restricted 
to time delays which are explicitly related to spatial 
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PILOT 

SIGNAL 



Figure 8. Adaptive Signal-Matching System 
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O o 



direction. These assumptions are not overly restrictive 
in normal applications; however, they may be eliminated by 
application of a more general theory of space-time cross- 
correlation functions. Such a theory has been suggested by 
Childers and Reed [5] . Use of these concepts has been con- 
sidered by Cheng and Tseng [4] in the formulation of 
signal-to-noise ratio for arrays. 
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